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The main objective of this article are two-fold. First, we introduce some general principles on phase 
transition dynamics, including a new dynamic transition classification scheme, and a Ginzburg- 
Landau theory for modeling equilibrium phase transitions. Second, apply the general principles 
and the recently developed dynamic transition theory to study dynamic phase transitions of PVT 
systems. In particular, we establish a new time-dependent Ginzburg-Landau model, whose dynamic 
transition analysis is carried out. It is worth pointing out that the new dynamic transition theory, 
along with the dynamic classification scheme and new time-dependent Ginzburg Landau models for 
equilibrium phase transitions can be used in other phase transition problems, including e.g. the 
ferromagnetism and superfluidity, which will be reported elsewhere. In addition, the analysis for 
the PVT system in this article leads to a few physical predications, which are otherwise unclear 
from the physical point of view. 

PACS numbers: 



I. INTRODUCTION 

The main objectives of this article are 

1) to introduce a new dynamic phase transition clas- 
sification scheme for nonlinear systems, 

2) to present a new time-dependent Ginzburg-Landau 
model for general equilibrium phase transitions, 
and 

3) to study the dynamic phase transitions for PVT 
systems. 

As we know, there are several classification schemes for 
phase transitions. The one used most often is the classi- 
cal Ehrenfest classification scheme, which is based on the 
order of differentiability of the energy functional. The 
new dynamic phase transition classification scheme we 
propose here is based on the dynamic transition theory 
developed recently by the authors [HE], which has been 
used to study numerous problems in sciences and engi- 
neering, including in particular problems in classical and 
geophysical fluid dynamics, biology and chemistry, and 
phase transitions. The new scheme introduced in this 
article classifies phase transitions into three categories: 
Type-I, Type-II and Type-Ill, corresponding respectively 
to the continuous, the jump and mixed transitions in the 
mathematical dynamic transition theory. 

We shall see from the applications to PVT systems 
that this new classification scheme is more transparent. 
In addition, it can be used to identify high-order (in the 
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Ehrenfest sense) transitions, which are usually hard to 
derive both theoretically and experimentally. 

Second, based on the le Chatelier principle and some 
general characteristics of pseudo-gradient systems, we in- 
troduce a general principle, leading to a unified approach 
to derive Ginzburg-Landau type of time-dependent mod- 
els for equilibrium phase transitions. In view of this new 
dynamic modeling and the dynamic transition theory, we 
shall see that the states after transition often include not 
only equilibria, but also some transient states, which are 
physically important as well. 

Third, with the aforementioned principles in our dis- 
posal, we derive a new dynamic model for PVT system, 
whose dynamic transition can be explicitly obtained, 
leading to some specific physical predictions. 

As we know, a PVT system is a system composed 
of one type of molecules, and the interaction between 
molecules is governed by the van der Waals law. The 
molecules generally have a repulsive core and a short- 
range attraction region outside the core. Such systems 
have a number of phases: gas, liquid and solid, and a solid 
can appear in a few phases. The most typical example of 
a PVT system is water. In general, the phase transitions 
of PVT systems mainly refer to the gas-liquid, gas-solid 
and liquid-solid phase transitions. These transitions are 
all first order in the Ehrenfest sense (i.e., discontinuous) 
and are accompanied by a latent heat and a change in 
density. 

A PT-phase diagram of a typical PVT system is 
schematically illustrated by Figure [l] where point A is 
the triple point at which the gas, liquid, and solid phases 
can coexist. Point C is the Andrews critical point at 
which the gas-liquid coexistence curve terminates [3]. 
The fact that the gas-liquid coexistence curve has a crit- 
ical point means that we can go continuously from a 
gaseous state to a liquid state without ever meeting an 
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observable phase transition, if we choose the right path. 
The critical point in the liquid-solid coexistence curve has 
never been observed. It implies that we must go through 
a first-order phase transition in going from the liquid to 
the solid state. 




In addition, in the gas-solid and liquid-solid transi- 
tions, there also exist metastable states. For the gas-solid 
case, the metastable states correspond to the superheated 
solid and supercooled liquid, and for the liquid-solid case, 
the metastable states correspond to the superheated solid 
and supercooled liquid; see Physical Conclusion V.2. 

This article is organized as follows. The new dynamic 
classification scheme and the Ginzburg-Landau model for 
equilibrium phase transitions are given in Section II. Sec- 
tion III introduces a time-dependent model for PVT sys- 
tems, which is analyzed in Section IV. Section V provides 
physical consequences of the analysis and some new phys- 
ical predictions. The abstract dynamical transition theory 
is recapitulated in the Appendix A for convenience. 



FIG. 1: Coexistence curves of a typical PVT system: A is 
the triple point, C is the critical point, and the dashed cure 
is a melting line with negative slope. 

The phase diagram in the PV-plane is given by Figure 
[2] where the dashed curves represent lines of constant 
temperature. In the region of coexistence of phases, the 
isotherms (dashed lines) are always flat, indicating that 
in these regions the change in volume (density) occurs 
for constant pressure p and temperature T. 
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FIG. 2: PF-phase diagram, the dashed line represent 
isotherms. 

One of the main objectives of this article is to study 
the dynamic phase transitions beyond the Andrews point 
C, and to make some physical predications. 

For this purpose, first, by both mathematical and phys- 
ical considerations, a proper form of the Gibbs free en- 
ergy functional is given, and then the time-dependent 
Ginzburg-Landau type of model for PVT systems is de- 
rived naturally using the general principle as mentioned 
above. 

Then, using the dynamical transition theory, we derive 
a system of critical parameter equations to determine the 
Andrews point C = [Tq,Pc) given by (V.6l. Then we 
show that the gas-liquid transition can be either first- 
order, second order or third order, as demonstrated in 
Physical Conclusion V.l, if we go from a gas to a liquid by 
a path without passing through the gas-liquid coexistence 
curve. 



II. GENERAL PRINCIPLES OF PHASE 
TRANSITION DYNAMICS 

In this section, we introduce a new phase dynamic 
transition classification scheme to classify phase transi- 
tions into three categories: Type-I, Type-II and Type-Ill, 
corresponding mathematically continuous, jump mixed 
transitions, respectively. 

Then a new time-dependent Ginzburg-Landau theory 
for modeling equilibrium phase transitions in statistical 
physics is derived based on the le Chatelier principle and 
some mathematical insights on pseudo-gradient systems. 



A. Dynamic Classification Scheme 

In sciences, nonlinear dissipative systems are gener- 
ally governed by differential equations, which can be ex- 
pressed in the following abstract form 



du 
~dt ~ 
«(0) 



L\u + G(u, A), 
= <P, 



(11.1) 



where u : [0, 00) — > X is the unknown function, A S 
(JV > 1) is the control parameter, X and X± are two 
Banach spaces with X\ C X being a dense and compact 
inclusion, L\ = —A + B\ and G(-,A) : X± — > X are 
C r (r > 0) mappings depending continuously on A, L\ : 
X\ — > X is a sectorial operator, and 



A : X\ —> X a linear homeomorphism, 
B\ : X\ — > X a linear compact operator. 



(II.2) 



In following, we introduce some basic and universal 
concepts in nonlinear sciences. 



First, a state of the system (II.l I at A is usually referred 
to as a compact invariant set In many applications, 
£a is a singular point or a periodic orbit. A state £a 
of (II.l I is stable if £a is an attractor, otherwise £a is 



called unstable. 

Second, we say that the system (II.l I has a phase tran- 
sition from a state Sa at A = Aq if Sa is stable on A < Ao 
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(or on A > Ao) and is unstable on A > Ao (or on A < A ). 
The critical parameter Ao is called a critical point. In 
other words, the phase transition corresponds to an ex- 
change of stable states. 

The concept of phase transition originates from the 
statistical physics and thermodynamics. In physics and 
chemistry, " phase" means the homogeneous part in a het- 
erogeneous system. However, here the so called phase 
means the stable state in the systems of nonlinear sci- 
ences including physics, chemistry, biology, ecology, eco- 
nomics, fluid dynamics and geophysical fluid dynamics, 
etc. Hence, here the content of phase transition has been 
endowed with more general significance. In fact, the 
phase transition dynamics introduced here can be ap- 
plied to a wide variety of topics involving the universal 
critical phenomena of state changes in nature in a unified 
mathematical viewpoint and manner. 



Third, if the system (II.ll possesses the gradient- type 
structure, then the phase transitions are called equilib- 
rium phase transition; otherwise they are called the non- 
equilibrium phase transitions. 

Fourth, classically, there are several ways to classify 
phase transitions. The one most used is the Ehren- 
fest classification scheme, which groups phase transitions 
based on the degree of non-analyticity involved. First or- 
der phase transitions arc also called discontinuous, and 
higher order phase transitions (n > 1) are called contin- 
uous. 

Here we introduce the following notion of dynamic clas- 
sification scheme: 



Defi nitio n II. 1 Let Ao G M. N be a critical point of (II. 1), 
and (II. 1) undergo a transition from state Yi\ toT? x . 
There are three types of phase transitions for |7J. 1 ) at 
A = Ao, depending on their dynamic properties: contin- 
uous, jump, and mixed as given in Theorem \A.1\ which 
are called Type-I, Type-II and Type-Ill respectively. 

The main characteristics of Type-II phase transitions 
is that there is a gap between and E| at the critical 
point Ao- In thermodynamics, the metastable states cor- 
respond in general to the super-heated or super-cooled 
states, which have been found in many physical phenom- 
ena. In particular, Type-II phase transitions are always 
accompanied with the latent heat to occur. 

In a Type-I phase transition, two states Yi\ and S| 
meet at Ao, i.e., the system undergoes a continuous tran- 
sition from to S^. 

In a Type-Ill phase transition, there are at least two 
different stable states and £3 at Ao, and system un- 
dergoes a continuous transition to or a jump transi- 
tion to S3, depending on the fluctuations. 

It is clear that a Type-II phase transition of gradient- 
type systems must be discontinuous or the zero order 
because there is a gap between T,\ and T, 3 X (2 < j < K). 
For a Type-I phase transition, the energy is continuous, 
and consequently, it is an n-th order transition in the 
Ehrcnfest sense for some n > 2. A Type-III phase tran- 
sition is indefinite, for the transition from to Yi\ it 



may be continuous, i.e., dF jd\ — dF^~ /dX, and for 
the transitions from to S| it may be discontinuous: 
dF- /d\ ^ dF+ /d\ at A = A . 

Finally, we recall the physical meaning of the deriva- 
tives of F(u, A) on A in statistical physics. In thermody- 
namics, the parameter A generally stands for temperature 
T and pressure p 

X = (T,p). 

The energy functional F(u, T, p) is the thermodynamic 
potential, and u is the order parameter. The first order 
derivatives are 



OF 
df 



= S 



dF 
dp 



= V 



the entropy, 

the phase volume, 



and the second order derivatives are 



1 d 2 F 



T-^p^ — C p the heat capacity in constant pressure, 



k the compression coefficient, 



V dp 2 
1 d 2 F 
VdTdp 



— a the thermal expansion coefficient, 



for the PVT system. 

Thus, for the first order phase transition, the disconti- 
nuity of dF/dp at critical point Ao = (To,po) implies the 
discontinuity of phase volume: 



a jrl 8F+ dF- 



AV = V z - V 1 = 



dp dp ' 



(II.3) 



and the discontinuity of dF/dT implies that there is a 
gap between both phase entropies: 



AS = S 2 - S 1 = - 



dF+ dF~ 



dT dT 



(II.4) 



Physically, AS is a non-measurable quantity, hence, we 
always use latent heat AH = TAS to determine the first 
order phase transition, and AH stand for the absorbing 
heat for transition from phase T,\ to phase In the 
two phase coexistence situation, AH and AV are related 
by the Clausius-Clapeyron equation 



dp 
dT 



AH 
TAV 



For second order phase transitions the variance of heat 
capacity (or specific heat), compression coefficient, and 
the thermal expansion coefficient at critical value are im- 
portant static properties in equilibrium phase transition, 
which are measurable physical quantities. 
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B. New Ginzburg-Landau Models for Equilibrium 
Phase Transitions 



In this subsection, we introduce the time-dependent 
Ginzburg-Landau model for equilibrium phase transi- 
tions. 

We start with thermodynamic potentials and the 
Ginzburg-Landau free energy. As we know, four ther- 
modynamic potentials- internal energy, the enthalpy, the 
Hclmholtz free energy and the Gibbs free energy-are use- 
ful in the chemical thermodynamics of reactions and non- 
cyclic processes. 

Consider a thermal system, its order parameter u 
changes in Q, C R" (1 < n < 3). In this situation, the 
free energy of this system is of the form 



must be the minimal value point of the thermodynamic 
potential. 

By the mathematical characterization of gradient 
systems and the le Chatelier principle, for a system 
with thermodynamic potential TL(u 1 A), the governing 
equations arc essentially determined by the functional 
7i(u,\). When the order parameters (u\,--' ,u m ) are 
nonconserved variables, i.e., the integers 



Ui(x,t)dx = a,i(i) ^ constant. 



then the time-dependent equations are given by 



H(u,\)=H 



. m 



Hi\Viii\ 2 + g(u,Vu,X) 



dx 



8U,, 

~di 
du 
dn 

u(x,0) = if(x) 



\dn 



fa— H(u,\) +$i(u,Vu,A), 
(oru\ dn = 0), 



(11.7) 



(II.5) 

where N > 3 is an integer, u — (iti,--- ,u m ),Ui = 
/i;(A) > 0, and g(u, Vit, A) is a C r (r > 2) function of 
(it, Vm) with the Taylor expansion 



for any 1 < i < m, where 5 /Sui are the variational deriva- 
tive, Pi > and <&i satisfy 



Jn : dm 



x = 0. 



(II. 



N 

g(u, Vu, A) = V a ljk u t D jUk + V a^u 1 + o(\u\ N )- fX, The condition (|IL8) is required by the Le Chatelier prin 

|7|=1 



jn.6) 

where / = (ii, • • • ,i m ),ik > are integer, |/| = J2T=i 
the coefficients and a/ continuously depend on A, 
which are determined by the concrete physical problem, 
u 1 = Ui ■ ■ ■ u % ™ and fX the generalized work. 

Thus, the study of thermal equilibrium phase transi- 
tion for the static situation is referred to the steady state 
bifurcation of the system of elliptic equations 



ciple. In the concrete problem, the terms $i can be de- 
termined by physical laws and (II 



When the order parameters are the number density 
and the system has no material exchange with the exter- 
nal, then Uj (1 < j < m) are conserved, i.e., 



Uj(x, t)dx — constant. 



(II.9) 



5u 

du 

On 



This conservation law requires a continuous equation 
du 



H(u,X) = 0, 

0, (or u\ d n 



dt 



= -V- Jj{u,X), 



(11.10) 



0), 



where Jj(u, A) is the flux of component Uj. In addition, 
Jj satisfy 



where S/Su is the variational derivative. 

A thermal system is controled by some parameter A. 
When A is for from the critical point Ao the system lies on 
a stable equilibrium state Si, and when A reaches or ex- 
ceeds Ao the state Si becomes unstable, and meanwhile 
the system will undergo a transition from Si to another 
stable state S2. The basic principle is that there often 
exists fluctuations in the system leading to a deviation 
from the equilibrium states, and the phase transition pro- 
cess is a dynamical behavior, which should be described 
by a time-dependent equation. 

To derive a general time-dependent model, first we re- 
call that the classical le Chatelier principle amounts to 
saying that for a stable equilibrium state of a system S, 
when the system deviates from S by a small perturbation 
or fluctuation, there will be a resuming force to retore 
this system to return to the stable state S. Second, we 
know that a stable equilibrium state of a thermal system 



(11.11) 



where /i; is the chemical potential of component ui, 
5 



E 



Mi = 



5u 



-H(u,\) - (j>j{u, Vw,A), (11.12) 



and (j}j(u, A) is a function depending on the other compo- 
nents Ui (i ^ j). When m = 1, i.e., the system consists 
of two components A and B, this term <pj = 0. Thus, 
from (11.10 )-( 11.12 1 we obtain the dynamical equations 



as follows 

duj_ 
dt 
du 

on on 
u(x,0) = (p(x), 



P 3 A 
Ian = 0, 



7i(u, A) — (j)j(u, V«, A) 



dAu 



Ian = 0, 



(11.13) 
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for 1 < j < m, where /3, > are constants, (/>,- satisfy 



Jn „■ 



6 u 



-H(u, \)dx = 0. 



(11.14) 



If the order parameters (u\ , • • • , Uk) are coupled to the 
conserved variables (u/ £+1 , • • • ,u m ), then the dynamical 
equations are 

-P iT7 -H(u,X)+$i(u,Vu,X), 



dt 

dt 
duj 
dn 



8u 



H(u, A) — 4>j (u, Vm, A) 



|an = (orui|ao = 0), 



(11.15) 



du 3 1 n 
-5— \an = u, 
on 

u(x, 0) = ip(x). 
for 1 < i < k and k - 



dAuj 
dn 



\an = 0, 



1 < j < m. 



The model (11.15 1 gives a general form of the governing 



equations to thermodynamic phase transitions. Hence, 
the dynamics of equilibrium phase transition in statistic 
physics is based on the new Ginzburg-Landau formula- 
tion (III. 15 I. 



Physically, the initial value condition u(0) — ip in 
(11.15 1 stands for the fluctuation of system or pertur- 



bation from the external. Hence, ip is generally small. 
However, we can not exclude the possibility of a bigger 
noise (p. 



From conditions (II. 8 1 and (11.14) it follows that a 



(11.16) 



steady state solution uo of (11.15) satisfies 

$i(u , Vu ,X) = yi<i<k, 
A<j}j(uo, Vu , A) = VA; + l<j<m. 

Hence a stable equilibrium state must reach the minimal 
value of thermodynamic potential. In fact, uq fulfills 

/3i—H(u ,\)-<$>i(u ,\/u ,X) = 0, 



/JjA-j — H(u , A) — A(f>j(u , V«o, A) - 0, 



duj 
dn 



dn 



Su 



an 



\dn 



(orui| ao = 0), 



(11.17) 



dAu J 
dn 



0, 



for 1 < i < k and k + 1 < j < m. Multiplying 
<I>j(ito, Vito, A) and 7 (mq, Vuq, A) on the first and the 
second equations of ( |II.17 l respectively, and integrating 
them, then we infer from (II. 8 1 and ( |II.14 | that 



/ Y^^(u O: \7u Ol \)dx = 0, 
Jn -. 

I V|V^(u ,Vu ,A)| 2 d2; = 0, 
Jn , 



which imply that (11.16 1 holds true 



III. TIME-DEPENDENT MODEL FOR PVT 
SYSTEMS 

In this section, we use the general principles derived 
in the last section to derive a time-dependent Ginzburg- 
Landau model for PVT systems, which will be used to 
carry out dynamic transition analysis and physical pred- 
ications in the next section. 



A. van der Waals equations and Gibbs energy 

The classical and the simplest equation of state which 
can exhibit many of the essential features of the gas- liquid 
phase transition is the van der Waals equation. It reads 



RT\ 



P J 



\v z + -v = 0, 



P 



P 



(III.l) 



where v is the molar volume, p is the pressure, T is the 
temperature, R is the universal gas constant, 6 is the re- 
vised constant of inherent volume, and a is the revised 
constant of attractive force between molecules. If we 



adopt the molar density p = 1/v to replace v in (III.l I, 
then the van der Waals equation becomes 



(bp + RT)p + ap 2 - abp 3 +p = 0. 



(III.2) 



Now, we shall apply thermodynamic potentials to in- 
vestigate the phase transitions of PVT systems, and we 
shall see later that the van der Waals equation can be 
derived as a Euler-Langrange equation for the minimiz- 
ers of the Gibbs free energy for PVT systems at gaseous 
states. 

Consider an isothermal-isopiestic process. The ther- 
modynamic potential is taken to be the Gibbs free en- 
ergy. In this case, the order parameters are the molar 
density p and the entropy S, and the control parameters 
are the pressure p and temperature T. The general form 
of the Gibbs free energy for PVT systems is given as 



G(p, S, T, p) 



^ Ul \Vp\ 2 + 1 -u 2 \VS\ 2 + g(p,S,T,p) 



ST — a(p,T,p)p dx 



(111.3) 



where g and a are differentiable with respect to p and 
S, O C R 3 is the container, and ap is the mechanical 
coupling term in the Gibbs free energy, which can be 
expressed by 



a(p,T,p)p= pp- \bp 2 p 



(III.4) 



where 6 = b(T, p) depends continuously on T and p. 

Based on both the physical and mathematical consid- 
erations, we take the Taylor expansion of g(p, S,T,p) on 
p and S as follows 



1 



-aip 



0k 



P2 Sp 2 



1 



1 



a 2 p i + -azp\ (III.5) 
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where ctj (1 < i < 3), /3i and /?2 depend continuously on 
T and p, and 



« t = cti(T,p) > i = 2,3, 
j 9i = /3i(T,p)>0. 



(III.6) 



B. Dynamical equations for PVT systems 

In a PVT system, the order parameter is u = (p, S), 

P = Px-po, S = Si-S , 

where pi and Si(i = 0,1) represent the density and en- 
tropy po,So are reference points. Hence the conjugate 
variables of p and S are the pressure p and the temper- 
ature T. Thus, by the standard model (II.7), we derive 



from (III. 3 )-( III.5 1 the following general form of the time- 



dependent equations governing a PVT system: 
dp 

at 



PiAp - (ai + bp)p + a 2 p 2 



95 
~dt 



- a 3 p 3 - 2f3 2 pS + P , 
= fj, 2 AS - piS - fop 2 + T. 



(HI.7) 



Although the domain f2 depends on T and p, we can 
still take the Neumann boundary condition 



dp 
dn 



<>{> 



= 0, 



as 

dn 



= 0. 



(III. 



Oil 



An important special case for PVT systems is that the 
pressure and temperature functions are homogeneous in 
O. Thus we can assume that p and S are independent of 
x £ ft, and the free energy pll.3| with ( |III.4[ ) and ( |III.5[ ) 
can be expressed as 



G(p,S,T,p) =^«iP 2 
1 



+ \f3iS 2 + I3 2 S P 2 



1 



;Ci 2 p 



1 



-,u 3 p 



bp 2 p - pp- ST. (111.9) 



From (III. 9 1 we get the dynamical equations as 
dp 



dt 

dS^ 
~dt 



= —(ai + bp)p + a 2 p 
a 3 p 3 - 2/3 2 Sp + p, 
-f3 1 S-f3 2 p 2 + T. 



(111.10) 



Because (3i > for all T and p, we can replace the second 
equation of (III. 10) by 

^\T-p 2 p 2 ). (III.ll) 



S 



dp 
dt 



Then, (|III.10|) are equivalent to the following equation 

(a 1 + bp + 2f3i 1 (3 2 T)p+a 2 p 2 

-(a 3 -2p 2 l3^)p 3 +p. (111.12) 

It is clear that if a,\ — 0,2(3^ 1 p 2 — R,a 2 — a , (ax — 
2/3 2 2 f3^ 1 ) = ab, then the steady state equation of (III. 12) 
is referred to the van der Waals equation. 



IV. PHASE TRANSITION DYNAMICS FOR 
PVT SYSTEMS 



In this section we use (III. 12 1 to discuss dynamical 



properties of transitions for PVT systems, and remark 
that similar results can also derived using the more gen- 



eral form of the time-dependent model (|HL7 1 

Let po be a steady state solution of 
the transformation 

P = Po + p' ■ 



11.121. We take 



Then equation (III. 12) becomes (drop the prime) 
dp 



dt 



= Xp + a 2 p 



d3P 



(IV.l) 



where 



A = 2a 2 p Q - 3a 3 p 
a 2 = a 2 - 3a 3 p , 
a 3 = a 3 -2f3 2 Pi 1 . 



a 1 -bp-2fo^ 1 T, 



In the PT-plane, near a non-triple point A = (T*,p*) 
(see Figure [T]) , the critical parameter equation 



A = A(I»=0 



\T-T*\<8, \p-p*\<6 



for some S > 0, defines a continuous function T = <p[p), 
such that 



<0 if T>0(p), 
= ifT = 0(p), 
>0 if T<cf)(p). 



(IV.2) 



The main result in this section is the following dynamic 
transition theorem for PVT systems. 

Theorem IV.l Let Tq = <fi(po) and a 3 > 0. Then the 
system (IV.l) has a transition at (T,p) — (To,po), and 
the following assertions hold true: 



1. If the coefficient a 2 = a 2 (T,p) in {IV.l) is zero at 
(Tq,Po), i.e., a 2 (T ,p ) = 0, then the transition is 
of Type-I, as schematically shown in Figure [?] (a) 
and (b). 

2. If o, 2 (Tq,pq) ^ 0, then the transition is of Type-Ill 
(i.e., the mixed type), and the following assertions 
hold true: 

(a) There are two transition solutions near 
(T , p ) as 



p ± (T,p) = JL (a 2 ± v /ai + 4a 3 A 



(IV.3) 



(b) There is a saddle-node bifurcation at (T\,pi), 
where T\ > T and pi < p if <f>' (pv) > 0, and 
Po > Pi if <P'(Po) < 0. 
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(c) For <fr'(po) > 0, when a 2 (T ,po) > the 
transition diagrams are illustrated by Figure ^ 
(a)-(b), where p + is stable for all (T,p) near 
{To, Pa), and p — is stable, p~ a saddle for 
To < T < T\ and p\ < p < po, and p~ is 
stable, p = a saddle for T < To,p > po; 

(d) When a2(To,po) < 0, the transition diagrams 
are illustrated by Figure^a)-(b), where p~ is 
stable for all (T,p) near (Tq,pq), and p = 
is stable, p + a saddle for Tq < T < T± , p± < 
p < Po, and p + is stable, p = a saddle for 
T < T ,po < p. 





(b) 



FIG. 4: Type-Ill (mixed) transition for a2 > 0. 




T 



FIG. 3: Continuous transition for the case where <f>'(po) > 



Proof. By (IV.l I, the theorem follows from the tran- 



sition theorems (Theorems A. 2 and A. 3 1 and the singu- 
larity separation theorem (Theorem A.4|. We omit the 
detailed routine analysis. ■ 



PHYSICAL CONCLUSIONS AND 
PREDICTIONS 



To discuss the physical significance of Theorem |IV.l 



we recall the classical p^-phase diagram given by Figure 
[5] We take p — l/V to replace volume V, then the gas- 
liquid coexistence curve in the p-P plane is illustrated by 
Figure [6] 

According to the physical experiments, in the gas- 
liquid coexistence region in the p~p phase diagram, there 
exist metastable states. Mathematically speaking, the 
metastable states are the attractors which have a small 
basin of attraction. In Figure [6] the dashed lines aa! 
and bb' represent the metastable states, and the points 




FIG. 5: Type-Ill (mixed) transition for 02 < 0. 



in aa' correspond to super-heated liquid, the points in 
bb' correspond to super-cooled gas. The p — p phase 
diagram shows that along the isothermal line T > Tc 
where C = (Tc,Pc) is the Andrews critical point, when 
the pressure p increases, the density p varies continu- 
ously from gaseous to liquid states. However along the 
isothermal line T < Tc when the pressure p increases to 
Po the density p will undergo an abrupt change, and a 
transition from gaseous state a to a liquid state b accom- 
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liquid 



T<T C 




T>T C 



FIG. 6: The gas-liquid coexistence curve in p — p plane: point 
C is the Andrews critical point, the dashed lines ad and bb' 
stand for the metastable states. 



panicd with an isothermal exothermal process to occur. 
Likewise, such processes also occur in the gas-solid and 
liquid-solid transitions. 

We now return to discuss Theorem 
state solution p — p (T,p) of (III. 12 1 



WJ\ The steady 



can be taken to 
represent a desired state in investigating different transi- 
tion situation. For example, for studying the gas-liquid 
transition we can take po as the gas density near the tran- 
sition temperature and pressure, and for the liquid-solid 
transition we take po as the liquid density, i.e., the lower 
state density. For convenience, in the following we alway 
consider the gas-liquid transition, and take po as the gas 

density. 

Let T = 0Oq), and <j>'(p ) > 0. Then Theorem |IVT 
implies the following physical conclusions: 

First, near (To,po), there are three stable equilibrium 
states of (III. 12 1 for each of the three cases: a 2 > 0, 



a 2 = 0, or a 2 < 0: 

a) If a 2 < 0, they are 

ip + = Po + P + for \T-T \ <Sk \p- Po \ <5, 
<p° = po for T < T or p < p , 

= po + p~ for T < T & p > p ; 



b) If a 2 = 0, they are 

<^ ± = Po + P ± for T < T & p > po, 
tp° = po forT > T orp< p , 



(V.l) 



(V.2) 



where p ± = 
c) If a 2 < 0, they are 

ip + = po + p + for T < T & p > po, 

V a = po for T > T or p < p , (V.3) 

■p-= P o + p- for |r-T | <5k |p — po| <*■ 

Here the state ip + represents the liquid density, ip° the 
real gas density, and ip~ the underlying state density. 



Second, putting S = j3^(T - (3 2 p 2 ) into ( |lll.9| , the 
free energy becomes 

G(p,T,p) = - l - ( 3- 1 T 2 + l -{a 1 + bp+2fc 1 p 2 T)p 2 

- l -a 2 p*+ -^a 3 -2Plp- 1 )p A -pp. 



The values of G at equilibrium states of (IV.l I are given 

by 

G( P o + ft =G{ P o) - |(p ± ) 3 + f (p^f 



=G(p )~(p ± ) 2 -|(/> ± ) 3 . 



(V.4) 



By A = at (T Q ,p ), it follows from |vl])-fO| that for 
(T,p) near (To,po) we have 

G(ip+) < G(<p°) < G(ip-) for a 2 > 0, 

G(ip+) = G(ip~) < G(tp°) for a 2 = 0, (V.5) 

G{tp~) < G(cp°) < G(cp+) for a 2 < 0. 



Third, according to Theorem |I V . 1 1 and the transition 
diagrams (Figures 3 51, there are two types of transition 
behaviors characterized by some functions $ = p(T,p) 
near (T, p) = (T ,po)i which are called the transition 
functions. Here, for simplicity we fixed p = po and con- 
sider $ = p(T) as function of T. 

If a 2 > 0, the transition functions are 



$+(T) = 



$"(T) 



Po(T) 



for T > T*, 



p (T) + p+(T) for T < T*, 
p (T) for T > T , 

Po(T)+p- forT<T , 



for some T <T* <T X . $+(T) has a finite jump at T = 
T* as shown in Figure [7^ a), and <i>~(T) is continuous, 
but has a discontinuous derivative at T = To; see Figure 

If a 2 =0, the functions are 



$+(T) 
$-(T) 



Po(T) for T > T , 

po{T)+p+(T) forT<T , 

Po(T) for T > T , 

p (T)+p-(T) forT<T , 



and $ + and are continuous with discontinuous 
derivatives at T = T ; see Figure [8ja) and (b). 
If a 2 < 0, the functions are 



p (T) for T > T , 

Po(T)+p+{T) forT<T , 

Po (T) for T>T*, 

p (T) + p-(T) for T<T\ 
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for some Tq < T * <T±, <f> + is continuous with a discon- 
tinuous derivative at T = Tq as shown in Figure |9ja) , 
and <I>~ is discontinuous at T = T* as shown in Figure 
§b). 

Fourth, based on physical facts, near the gas-liquid 
transition point (Tq,Po) the density p is a decreasing 
function of T, therefore the transition process <I> + is real- 
istic, and <1>~ is unrealistic. From the transition function 
$ + (T), the gas-liquid transition can be well understood. 



T* 




To 
(b) 



T 
(a) 




FIG. 8: The transition functions for 0,2 = 0: (a) the curve of 
p = $ + , (b) the curve of p = $~. 



P+P 



(a) 



FIG. 7: The transition functions for 0,2 > 0: (a) the curve of 
p = $ + , (b) the curve of p = $~. 



Fifth, we now explain the gas- liquid transition by 
using the transition function <& + (T). Let the two curves 

X = X(T,p)=0 and a 2 = a 2 (T,p) = (V.6) 



intersect at C — (Tc,pc) in PT-plane (see Figure 10 1, 
and the curve segment AB of A = is divided into two 
parts AC and CB by the point C such that 

a 2 (T,p)>0 for (T,p) e AC, 
a 2 (T,p)<0 for {T,p)eCB. 

We shall see that the point C = (Tc,pc) is the Andrews 
critical point, and the curve segment AC is the gas- liquid 
coexistence curve, as shown in Figure [l] 

In fact, when (T ,p ) € AC, a 2 (T ,p ) > 
transition function <i> + (T) has a jump at T 
Figure |6ja)): 



and the 
= T* (see 



Pa -> Pa 



0-2 

— > 0. 

03 



Hence the system undergoes a transition from a gaseous 
state to a liquid state with an abrupt change in density. 



T 

(b) 



FIG. 9: The transition functions for < 0: (a) the curve of 
p = $ + , (b) the curve of p = $~. 



On the other hand, by (V.5) there is an energy gap be 



tween the gaseous and liquid states: 

AE = G(po + p + ) — G(po) < atT = T*. 

This energy gap \AE\ stands for a latent heat, and 
AE < shows that the transition from a gaseous state 
to a liquid state is an isothermal exothermal process, and 
from a liquid state to gaseous state is an isothermal en- 
dothermal process. 

When (T , Po ) = {T c ,p c ),a 2 {T c ,Pc) = and the 
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FIG. 10: 



transition function <I> + (T) is continuous as shown in Fig- 
ure Ma). Near T — T c , 



<S>+{T) = p {T) + ^/a7 3 for T<T C . 
By the Landau mean field theory, we have 

A(T) = a(Tc — T) (a > a constant). 



(V.7) 



Thus we infer from (V.4| and (V.7) that 

„2 



G(<S>+(T)) = G(p Q ) 



a 
-la- 



(T - T c ) for T <T C 



The difference of the heat capacity at T — Tc is 

AC = -T C ^ (G($+(T C + 0)) - G( Po )) = ^-T c > 0. 

Namely the heat capacity has a finite jump at T = Tc, 
therefore the transition at T — Tc is of the second order. 

When (Tq,pq) € CB,a,2(T ,po) < and the transition 
function <!> + (T) is continuous as shown in Figure [9j a). 
Near T = T , 

$+(T) = po(T) + ^A(T)+ (|A|) T<T . 

ml 



Notice that (IV. 3 1 implies that 

+ A 
P = — 

a-2 



o(\X\). 



Then we deduce from (V.4l and (V.7 1 that if T > T 



G(<J> + (T)) = G(p (T)), 



and if T < To, 
G($ + (T)) = G( Po (T)) 



6a a 



;(To-T) 3 + o(|T -T| 3 ). 



Namely, the free energy G($ + (T)) is continuously dif- 
fcrcntiablc up to the second order at T = To, and 
the transition is of the third order. It implies that as 
(T),po) G GB the Type-I transition at (T),po) can not 



be observed by physical experiments. Therefore we can 
derive the following physical conclusion: 

Physic al Co nclusion V.l.The point C — (Tc,pc) sat- 
isfying (V.6) corresponds to the Andrews critical point, 



at which the gas-liquid transition is of the second or- 
der. In addition, when A(T),po) — 5 a 2(To,po) > 0. 
the gas-liquid transition at (T ,p ) is of the first or- 
der accompanied with a latent heat to occur, and when 
X(Tq,pq) = 0,a2(To,po) < 0, the transition at (To,po) is 
of the third order. 

Sixth, as X(T ,po) = and a 2 (T ,p ) > 0, the gas- 
liquid transition point (T*,p*) is in the range To < T* < 
Ti and p\ < p* < p ; see Figure [4j a)- (b). In fact, in the 
region To < T < T\ and p\ < p < po, the two stable 
states </?° = p (T,p) and ip + — p (T,p) + p + (T,p) are 
attractors, each possessing a small basin of attraction. 
Therefore they correspond to metastable states, and <p° 
can be considered as a super cooled gas, while ip + can be 
considered a super heated liquid. 

Finally, likewise, we can also discuss the gas-solid 
and liquid-solid transitions, and derive the following 
physical conclusion: 

Physical Conclusion V.2. In the gas-solid and liquid- 
solid transitions, there also exist metastable states. For 
the gas-solid case, the metastable states correspond to 
the superheated solid and supercooled liquid, and for the 
liquid-solid case, the metastable states correspond to the 
superheated solid and supercooled liquid. 



APPENDIX A: RECAPITULATION OF THE 
DYNAMIC TRANSITION THEORY OF 
NONLINEAR SYSTEMS 

In this appendix we recall some basic elements of the 
dynamic transition theory developed by the authors [1 
2 , which are used to carry out the dynamic transition 
analysis for the PVT systems in this article. 

Let X and X\ be two Banach spaces, and IjCla 
compact and dense inclusion. In this chapter, we always 
consider the following nonlinear evolution equations 



du 
Itt 
u(0) 



= L x u + G(u, A), 



(A.l) 



where u : [0, oo) — > X is unknown function, and A 6 K 1 
is the system parameter. 

Assume that L\ : X\ — * X is a parameterized linear 
completely continuous field depending contiguously on 
A 6 M 1 , which satisfies 



L\ = —A + B\ a sectorial operator, 

A : X\ — » X a linear homeomorphism, 

B\ : X\ — > X a linear compact operator. 



(A.2) 
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In this case, we can define the fractional order spaces X a 
for cr e i 1 . Then we also assume that G(-, A) : X a — > X 
is C r (r > 1) bounded mapping for some < a < 1, 
depending continuously on A £ K 1 , and 



G(u,A)=o(||u|| x J, VAe 



(A.3) 



Hereafter we always assume the conditions ( |A.2[ ) and 
(A.3 1, which represent that the system (A.l I has a dissi- 



pative structure. 

In the following we introduce the definition of transi- 
tions for (A.l I. 



Definition A.l We say that the system \A.ly has a 
transition of equilibrium from (u, A) = (0, Ao) on A > Ao 
(or A < Ao) if the following two conditions are satisfied: 

(1) when A < Ao (or A > Ao),u — is locally asymp- 
totically stable for (A.l); and 



(2) when A > Ao (or A < Ao), there exists a neighbor- 
hood U C X of u = independent of A, such that 
for any tp £ U \ T\ the solution u\(t,p) of (A.l) 
satisfies that 

]imsnp\\u\(t,<p)\\x > S(X) > 0, 

t — >-OC 

lim <5(A) > 0, 



where T\ is the stable manifold of u = 0, with codim 
T A > 1 in X for A > A (or A < A ). 



Obviously, the attractor bifurcation of ( |A.1[ ) is a type 
of transition. However, bifurcation and transition are 
two different, but related concepts. Definition |A.l| defines 
the transition of (A.l I from a stable equilibrium point to 



other states (not necessary equilibrium state) . In general, 
we can define transitions from one attractor to another 
as follows. 

Definition A. 2 Let Y,\ d X be an invariant set of 



(A.l). We say that (A.l) has a transition of states from 
(Xa , Ao) on A > Ao for A < Ao) if the following condi- 
tions are satisfied: 

(1) when A < Ao (or A > Ao),Ea is a local minimal 
attractor, and 

(2) when A > Ao (or A < Ao), there exists a neighbor- 
hood U C X of T.\ independent of A such that for 



any tp £ U \ (T\ U Xa), the solution u(t, tp) of (A.l) 
satisfies that 

limsupdist(u(f, tp), Xa) > 8(X) > 0, 

t — >oo 

lim 5(A) = S > 0, 



where T\ is the stable manifolds of Xa with codim 

r A >i. 



Let the eigenvalues (counting multiplicity) of L\ be 
given by 



{&(A) eC | j 



1,2,- 



} 



Assume that 



Re 0i (A) 



< 
= 
> 



if A < A , 
if A = A , 
if A > A , 



VI < i < 



(A.4) 



Re /3j(X ) < 



Vj>m+1. (A.5) 



The following theorem is a basic principle of transitions 
from equilibrium states, which provides sufficient condi- 
tions and a basic classification for transitions of nonlinear 
dissipative systems. This theorem is a direct consequence 
of the center manifold theorems and the stable manifold 
theorems; we omit the proof. 



Theorem A.l Let the conditions (A.4) and (A.5) hold 
true. Then, the system (A.l) must have a transition from 
(u, X) — (0,Ao), and there is a neighborhood U C X of 
u = such that the transition is one of the following three 
types: 

(1) Continuous Transition: there exists an open 
and dense set U\ C U such that for any tp £ U\, 
the solution u\(t,p) of \A.lty satisfies 



lim limsup ||u,\(£, <p)\\x = 0. 



In particular, the attractor bifurcation of (A.l) at 
(0, Ao) is a continuous transition. 

(2) Jump Transition: for any A < A < A + e with 
some e > 0, there is an open and dense set U\ C U 
such that for any p G U\, 

limsup \\u\(t,tp)\\ x >S>0, 

t — >oo 

where S > is independent of X. This type of tran- 
sition is also called the discontinuous transition. 

(3) Mixed Transition: for any A < A < A + e with 
some e > 0, U can be decomposed into two open 
sets Ui and 11% (U^ not necessarily connected) : 



U = U£ 
such that 



U. 



2 > 



n U2 = 



lim limsup \\u(t, tp)\\x — "itp £ f/f 
limsup \\u(t, tp)\\ x > S > Wp £ C/ 2 A 

t — >oo 



The following theorem provides sufficient conditions 
for continuous transitions and gives local transition struc- 
ture. 
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u = 





(a) 



v. 



t 

(b) 



FIG. 11: Topological structure of the jump transition of ( A.l I 
when fc=odd and a > 0: (a) A < Ao; (b) A > Ao. Here the 
horizontal line represents the center manifold. 



FIG. 12: Topological structure of the continuous transition of 
(A.l) when fc=odd and a < 0: (a) A < Ao; (b) A > A . 



We consider the transition of (A.l ) from a simple crit- 
ical eigenvalue. Let the eigenvalues f3j(X) of L\ satisfy 



{< if A < A , 

= if A = A , 

> if A > A , 
Re/3j(A ) < 



(A.6) 



Vj > 2, 



where /?i(A) is a real eigenvalue. 

Let ei(A) and e*(A) be the eigenvectors of L\ and L* x 
respectively corresponding to /3i(A) with 

L\ ei = 0, L* Xo e{ = 0, < ei, e* >= 1. 



Let Q(x, A) be the center manifold function of (A.l 
A = Aq. We assume that 



near 



< G(xei + $(x, A ), Aq), e\ >= ax k + o(\x\ k ), (A.7) 

where k > 2 an integer and a/0a real number. 
We have the following transition theorems. 



Theorem A. 2 Under the conditions (A.6) and (A.l), if 



k=odd and a in (A.7) then the following assertions 
hold true: 



(1) If a > 0, then (A.l) has a jump transition from 
(0,Ao), and bifurcates on A < Ao to exactly two 
saddle points v x and v% with the Morse index one, 
as shown in Figure [77) 



(2) If a < 0, then (A.l) has a continuous transition 
from (0,Ao), which is an attractor bifurcation as 
shown in Figure {T2^ 

(3) The bifurcated singular points v± and v£ in the 
above cases can be expressed in the following form 



±|/3 1 (A)/a| 1 /'=-ie 1 (A)+ (|/3 i r/ fe - 1 ). 



Theorem A. 3 Under the conditions (A.6) and (A.7), if 



k=even and a ^ 0, then we have the following assertions: 



1. (A.l) has a mixed transition from (0, Ao). More 
precisely, there exists a neighborhood U C X of u = 
such that U is separated into two disjoint open 
sets Ui and C/^ by the stable manifold T\ of u = 
satisfying the following properties: 

(a) U = U? + C/ 2 A + T\, 

(b ) the transition in is jump, and 

(c) the transition in U£ is continuous. The local 



transition structure is as shown in Figure 13 



2. (A.l) bifurcates in to a unique singular point 
v A on A > Ao, which is an attractor such that for 
any (peU$, 

lim \\u(t,p)-v x \\ x = 0, 

t— >oo 



where u(t,ip) is the solution of (A.l) 



3. (A.l) bifurcates on A < Ao to a unique saddle point 
v A with the Morse index one. 

4- The bifurcated singular point v x can be expressed as 

v x = -(/MAj/a) 1 /**- 1 ^! + oflftf/C*- 1 )). 



We consider the equation (A.l I defined on the Hilbert 
spaces X = H,Xi = H±. Let L\ = —A + XB. For L\ 
and G(-,A) : H\ — » H, we assume that A : Hi — > H is 
symmetric, and 

< Au,u > H > c\\u\\ 2 Hi/2 , (A.8) 

< Bu, u > H > c\\u\\ 2 H , (A. 9) 

< Gu,u> H < - Cl |M|^ + c 2 |H|^, (A.10) 

where p > 2, c, c\, ci > are constants. 
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v'- 


u = 



Theorem A. 4 Assume the conditions (A. 3), \ A.J^ and 
(A.8)-(A.10), then (A.l) has a transition at (u, A) = 



(0, A ), and the following assertions hold true: 



(1) If u = is an even-order nondegenerate singular 
point ofL\ + G at X = Xq, then (A.l) has a singular 
separation of singular points at some (ui, Ai) £ Hx 
(-oo, A ). 



(2) If m = 1 and G satisfies (A.l) with a > if k=odd 
and a ^ if k=even, then (A.l) has a saddle- 
node bifurcation 
with Ai < Aq. 



-even, then 
at some singular point (ui,Ai) 



FIG. 13: Topological structure of the mixing transition of 
(A.l I when fc=even and a / 0: (a) A < Aq; (b) A = Ao; (c) 



A > Ao. Here U\ is the unstable domain, 
domain. 
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